Impact of spin-orbit coupling on the Holstein polaron 
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We utilize an exact variational numerical procedure to calculate the ground state properties of a 
polaron in the presence of a Rashba-like spin orbit interaction. Our results corroborate with previous 
work performed with the Momentum Average approximation and with weak coupling perturbation 
theory. We find that spin orbit coupling increases the effective mass in the regime with weak electron 
phonon coupling, and decreases the effective mass in the intermediate and strong electron phonon 
coupling regime. Analytical strong coupling perturbation theory results confirm our numerical 
results in the small polaron regime. A large amount of spin orbit coupling can lead to a significant 
lowering of the polaron effective mass. 



I. INTRODUCTION 

In much of condensed matter (magnetism excepted), 
the spin and orbital components of an electron are 
treated as independent degrees of freedom. Nonetheless, 
the non-relativistic approximation to the Dirac equation 
leads directly to the so-called Thomas term in the effec- 
tive Hamiltonian, which can be written as a spin-orbit 
coupling tcrmi This coupling can play a significant role 
in the electronic structure of semiconductors and metals, 
as documented, for example, in Ref. (0). More recently, 
interest has grown because of the burgeoning possibili- 
ties in the so-called field of spintronics, where the spin 
degree of freedom is specifically exploited for potential 
applications^ Control of spin will require coupling to the 
orbital motion, and hence spin-orbit coupling may play 
a critical role in understanding and exploiting various 
properties of such systems. 

Spin-orbit coupling, as described by Rashba^ is ex- 
pected to be prominent in two dimensional systems that 
lack inversion symmetry, including surface states. Many 
such systems have now been identified, among which are, 
for example, surface alloys, Li/W(110)£ Pb/Ag(lll)£l 
and Bi/Ag(lll)4 In all of these systems the possibility 
of other interactions remains; in particular recent workS. 
has focussed on the electron-phonon interaction, in the 
presence of Rashba spin-orbit interactions. In the first 
reference of Ref. (@), for example, the effective mass due 
to the electron-phonon interaction was shown, in weak 
coupling, to be enhanced by the spin-orbit interaction. 

More recently, attention has focussed on the proper- 
ties of a single electron interacting with oscillator de- 
grees of freedom^ in the presence of Rashba spin-orbit 
coupling.— These authors utilized the so-called momen- 
tum average (MA) approximation^ to examine the prop- 
erties of a single polaron also in the presence of spin-orbit 
coupling, but for a tight-binding model. They found that 
the effective mass generally decreases as a function of 
spin-orbit coupling, Vg; however, in the weak electron- 
phonon coupling limit, there is initially an increase in 
effective mass, in agreement with Cappelluti et al.— The 
primary purpose of this work is to present exact solutions 



to this problem, using Trugman's method , 13 > 14 along with 
some modified algorithms^ so that we can span the en- 
tire parameter regime. It turns out that the MA method 
is fairly accurate over the entire parameter range, except 
for low phonon frequency. 

We also develop a strong coupling expansion, based on 
the Lang-Firsov transformation^ following Ref. (fl8| ). 
As in the straightforward Holstein model, strong coupling 
describes fairly well the small polaron regime. Finally, 
the adiabatic limit of the Holstein model with Rashba 
spin-orbit coupling has been described recently in Ref. 
(fl9h . following Refs. ((201 ) and (f2lh for the simple Holstein 
model. In the strict adiabatic limit Grimaldi finds an 
intermediate state (large polaron) with the lowest energy, 
for coupling strengths just below that required for small 
polaron formation, in the pre sence of spin-orbit coupling 
(see Figs. 1 and 2 in Ref. (|19l )). Our search for this state 
will also be described in the present work. 



The paper is organized as follows. We first introduce 
the model of study; following Ref. (fIT[ ) it is the Holstein 
model with additional Rashba spin-orbit coupling, writ- 
ten for a tight-binding formulation. We note some of the 
features of the non-interacting (with respect to phonons) 
model. Unlike the continuum limit ^ there is not a sin- 
gularity at the bottom of the band; however, for weak 
spin-orbit coupling, a singularity remains very close by 
in energy, and causes a significant enhancement in the 
density of states at the bottom of the band. In Section 
III we present our numerical results, along with those 
from the strong coupling expansion and from the MA ap- 
proximation. As mentioned above, the exact numerical 
results confirm the conclusions from Ref. HO). Finally, 
we examine the low phonon frequency and intermediate 
electron-phonon and spin-orbit coupling regimes, where 
both perturbative and MA approaches are suspect. We 
are unable to rule out the presence of an intermediate 
phase completely, but find that its occurrence is unlikely, 
once quantum fluctuations are included. We close with 
a summary. 
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II. MODEL 

The standard formulation for spin-orbit interaction 
uses two different types of electronic band structure. 
The first is free electron-like, which results in parabolic 
bands,— and the second is tight-binding, which results in 
a periodic momentum dependence. While it is essentially 
always the case that the latter tends to the former for low 
electron fillings, this is not quite true when a Rashba- 
type spin-orbit interaction term is present. As shown in 
Ref.@, for example, the ground state for a single electron 
consists of a degenerate ring around the T— point. This 
results in an electronic density of states with a square- 
root singularity at the bottom of the band. For a tight- 
binding model, however, Covaci et alAi pointed out that 
this is not the case. We will adopt a tight-binding for- 
mulation here, and examine this difference more closely 
in the next sub-section. 

To study the single polaron with spin-orbit interaction 
we use a tight-binding Hamiltonian with Rashba-type 
spin-orbit interaction^ and a Holstein-typeiS electron- 
phonon interaction. In real space the Hamiltonian is: 

H =-t ( C L C J> + c ], a c i,«) 

<i,j>,a=\i 

+ V S ( ic la a x Pc i+M - ic \, a ofci+l:,tt + h - c -) 

-gu E c t c M( 8 . + fl ')+ w EE a ' a " W 

where c\ s (cj. s ) is the creation (annihilation) for an elec- 
tron at site i with spin index s, a\ (ai) is the creation (an- 
nihilation) operator for a phonon at site i, and er"' 3 , Oy^ 
designate the (a, /3) component of the usual Pauli ma- 
trices. The sum over i is over all sites in the lattice, 
whereas < i,j > means only nearest neighbor hopping 
is included. Here, as the notation already suggests, we 
confine ourselves to nearest neighbor hopping only. The 
energy scales are the hopping integral t, the strength of 
the Rashba spin-orbit interaction, Vs, the coupling of the 
electron to the oscillator degrees of freedom gu>E, and the 
Einstein phonon frequency, oje- In what follows we write 
all energy scales in terms of the hopping integral, t, which 
hereafter is set to unity. The ground-state properties of 
the Holstcin model in one and two dimensions near the 
adiabatic limit have recently been studied in Rcfs. (fl5l ) 
and lfl6h. Normally spin is not considered, since this 
ground state is degenerate with respect to spin. As the 
Rashba spin-orbit interaction is turned on, however, the 
two- fold degeneracy will be lifted. 



non-interacting model: ground state and 
effective mass 



To examine this model in detail, we use a 2 x 2 matrix 
to describe the spin sector, and begin by excluding the 



phonon part of the Hamiltonian. The remaining Hamil- 
tonian is diagonalized through Bloch states in momentum 
space, written as 

Hq = Y ekC l,a C Ka + Y ^k-CTa^cj^Ck^ (2) 

k,a k,a,/3 

where ej. = — 2t[cos(k x ) + cos(fc J/ )] and fik'C = 
2V's[sin(fcj ) )cr 2 ; — s'm(k x )a y ] (we set the lattice spacing a 
equal to unity). Diagonalizing this 2x2 matrix, we get 
two bands, which we name the upper and lower Rashba 
bands. The eigenvalues and eigenstates are given by 



H $> ± = e fc ,±tf±, 



(3) 



with eigenvalues 



£k,± = -2t[cos(k x ) + cos{k y )}±2V s \Jsm 2 (ky) + sin 2 (k x ), 

(4) 

and eigenvectors 
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t sm(ky) - ism(k x ) t 
Cfc t ± / .-,,, „ , , c kl 



sin (k y ) + sin (k x ) 



|0). (5) 



In contrast to the model with parabolic bands, this model 
has a four-fold degenerate ground state located at k x = 
ky = ± arctan(-^),— which can be seen clearly from a 
contour plot of the lower Rashba band in FigfTJ There are 
also four saddle points near the energy minimum points, 
which are located at k x = 0, k y = ± arctan(^-) and k y = 
0,k x = ±arctan(i^). As V s increases, the separation 
between minimum points and saddle points is enhanced 
(see below, in Figl^b)). The ground state energy for 
H is given by E = -4ty/l + Vg/{2t 2 ). Similarly, the 
effective mass along the diagonal is 
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m v/TTTfTPi 5 )' 



(6) 



where mo = l/(4t) is the bare mass in the absence of 
spin-orbit interaction, and mgo is the effective mass due 
solely to spin-orbit interaction. Note that the effective 
mass decreases due to spin-orbit interaction. Below we 
will turn on the electron phonon interaction, and the 
ground state energy (effective mass) will be further low- 
ered (raised) due to polaronic processes. 



B. non-interacting model: electron density of 
states 



The non interacting electron density of states (DOS) 
is defined for each band, as 



D s (E) = J2 6 (E-e ks ) 



(7) 



with s = ±1. 
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FIG. 1. (color online) Contour plots for lower Rashba band 
with Vs/t = 0,0.5,1.0,5.0. For Vs = 0, there is only one 
energy minimum point at k x — k y = 0. For Vs > 0, 
there are four energy minimum points located at k x = k y = 
± arctan(-j^). For nonzero Vs, there are also four saddle 
points near the energy minimum points, which are located at 
k x = 0, k y = ±arctan(-^-) and k y = 0, kx = ± axctan(^). 
As Vs increases, the separation between minimum points and 
saddle points is increased (see Fig^b)). 
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FIG. 2. (a) Non-interacting density of states D-(E) near the 
bottom of the band for Vs/t — 0,0.5, 1.0. In the inset the 
density of states in the whole band is shown for the same 
parameters. Note that the divergence at the bottom of the 
band has been shifted to higher value.— (b) The separation 
between energy minimum points and saddle points as a func- 
tion of spin orbit interaction Vs /t. 

In the main frame of Figl^a) we show the low energy 
DOS for various values of the spin orbit interaction Vs; 
note that this involves only D_ (E) as the upper Rashba 
band exists only at higher energies. Furthermore, infor- 
mation concerning the upper Rashba band can always be 
obtained through the symmetry 

D + {E) = D.{-E). (8) 

Figl2I a ) shows that a divergence introduced by the 
spin orbit interaction exists at higher energy^ and not 
at the bottom of the band, as occurs for a parabolic 
dispersion^ This shift is due to the separation of the 
energy minima from the saddle points in k-space, as 
shown in Figl^b). The saddle point energy is given by 
E s& a — — 2£(1 + y/l + (V s /t) 2 ), which is very close to 
the minimum energy Eq even for sizeable Vs/t, as is ev- 
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ident from the figure. This proximity of the divergence 
serves to elevate the value of the DOS at the bottom 
of the band. With no spin orbit coupling this value is 
D±(E = E = -At) = l/(47rt) (V s = 0). With spin or- 
bit coupling, however, an expansion around the minimum 
energy E Q = + V§ /2t 2 yields a DOS value 



D- (E = Eq) 



V2 1 

TT V S 



(9) 



Thus a discontinuity occurs as the spin orbit coupling is 
changed from zero — the DOS immediately has a diver- 
gence at the bottom of the band which, for any non-zero 
value of Vs, shifts to slightly higher energy. The inset 
shows D_{E) over a wider energy range. Further details 
are provided in Appendix A. 



III. GROUND STATE ENERGY AND 
EFFECTIVE MASS 

When the electron phonon interaction is turned on, 
the ground state energy (effective mass) will be lowered 
(increased) due to polaron effects. To study the polaron 
problem numerically, we adopt a variational method out- 
lined by Trugman and coworkers ; 13 ' 14 which could deter- 
mine polaron properties in the thermodynamic limit ac- 
curately. This method was recently developed by Alver- 
mann et ali£ and Li et al» to study the polaron problem 
near the adiabatic limit. In this paper we will adopt the 
numerical techniques described in Ref. [l5| . 



A. Strong coupling theory 

To investigate the strong coupling regime of the 
Rashba-Holstcin model for a single polaron, we use the 
Lang-Firso v 17 i 18 unitary transformation H = e s He~ s , 
where S = g^2 i a ni^ cr {a i — a\), and obtain 



H = H +T 



with 



and 



(10) 



(11) 



+iVsJ2( c t a °f*+y,e x l x i+y 

i 

-cl a afc w xjx. l+£ - h.c), (12) 



,A+B _ 



where Xj = exp{g(a,i — a])}. Using 
e" 4 e B e _1 / 2 [" 4,B l, the hopping part of the Hamiltonian be- 



1,(7,5 

i 

A^f^ApnHp^p^p^-h-c], (i3) 

where = cxp(±ga,;). The unperturbed bare Hamil- 
tonian, H provides the zeroth order energy for the po- 
laron, and is already diagonal for the single electron sec- 
tor. The eigenvalues are given by E n = nuiE — 5 2 we, 
where n is the total number of phonons. Clearly the 
ground state has n = 0, but remains 2N-fold degenerate, 
since the electron can occupy any one of the N sites and 
it can have either spin up or spin down. If we consider 
the hopping term T as a perturbation and apply degener- 
ate perturbation theory to the 2N-fold degenerate ground 
state, we need to diagonalize a 2N x 2N matrix. A sim- 
pler approach is to recognize that the momentum k is a 
good quantum number, and if we transform the original 
problem into k-space, we need only solve a 2 x 2 matrix 
which mixes the spin sectors; this results in essentially 
Eq. (|2J), but with an extra band narrowing factor e~ 9 . 
Thus we obtain the first order perturbation correction to 
the energy as 



E, 



k± 



(14) 



and the result is the familiar band narrowing factor that 
occurs when Vs = 0. 

The eigenstates from degenerate perturbation theory 
are now simply Bloch-like states, as found in the 
non- interacting theory, Eq. ([5]). Thus the degeneracy 
is broken, and a comparatively narrower band is formed 
with a minimum at a non-zero wave vector in the lower 
Rashba band, as found in the non-interacting case. To 
find the second order correction to the ground state en- 
ergy, we proceed as in Ref. [l8| . and find 



N 



E 



(2) 



E E 

riTOT^0,ni ,ri2,...— 0,1,.. .oo i ~ 1 



(ni,n 2 , —n N \ ph g (cea\ e iT\^k,-)ei ® \0) P h 

— riTOT^E 



(15) 



where utot is the total number of phonons and ^k,- is 
given in Eq. ([S])- With details shown in the appendix, 
we obtain 



E 



(2) 



_4e-^^i[ /( 2.g 2 )- / (. 9 2 )]- e -^ / ( ff 2 )^, 

UJE OJe 

(16) 

OO 

where f(x) = £ I£ « e */x[l + l/x + 2/x 2 + ...] (see 

n=l 

Appendix). In some of the ensuing discussion, we will 
use the constant A, familiar as the effective mass en- 
hancement from weak coupling perturbation theory for 
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the interacting electron gas. Here we use the definition^ 
A = 2g 2 WE^f, since l/(47rf) is the value of the non- 
interacting electron density of states for Vs — at the 
bottom of the band. Note that our definition of A dif- 



fers from that in Ref.jllj or Ref.|18f; both use the more 
conventional average density of states, l/(8t). Thus the 
ground state energy, excluding exponentially suppressed 
corrections, is 



E, 



GS 



-2nt\(l 



r + 



{2ntX) 



(17) 



and there is a correction of order 1/ A 2 compared to the 
zcroth order result. Corrections in the dispersion enter 
in strong coupling only with an exponential supression. 
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B. Weak coupling theory 

In the weak electron-phonon coupling regime, does 
spin-orbit coupling suppress or enhance the "polaron ef- 
fect" due to the electron-ion coupling? Weak coupling 
calculations with a parabolic electron dispersion^ showed 
an increase in the effective mass, for example, as the spin- 
orbit coupling was increased. Here we perform weak cou- 
pling perturbation theory, as described in Ref.@, with 
the same definitions, except that the tight binding dis- 
persion is used to describe the non-interacting electrons, 
as outlined in the previous section. A straightforward 
calculation yields the self energy to first order in A as 



S W eak(w + lS) = TrXtUJl 



1 



k,s=± 



uj + iS - lu e - Sh 



(18) 



The effective mass can be obtained by the derivative of 
the self energy 



m wcak 

mso 



1 



d_ 

duj 



S W eak(w + id)\ L 



(19) 



Near the adiabatic limit (lje 0), by expanding Ek.- 
around Eq, as shown in the appendix for the calculation 
of the DOS, we obtain 



TO wcak 

m S o 



1 



V2Xt 
~Vs~' 



(20) 



which shows a diverging effective mass as the spin-orbit 
coupling increases. In fact, there is a discontinuity for 
Vs = 0, as the result is simply T ^" , ° k = 1 + A/2, and 
mso -> m = l/2t, as given by Eq. ©. Eq. d^HJ) will 
have a limited domain of validity, however, as we will see 
below. 



C. Numerical Results 

In Fig. [31 we show the ground state energy and the ef- 
fective mass correction as a function of electron-phonon 
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FIG. 3. (color online) (a) Ground state energy difference 
Eos - E vs. A for V s /t = 0, 0.5, 1.0 and uj E /t = 1.0. Exact 
numerical results are compared with those from weak coupling 
perturbation theory (labeled "Pert." in the figure) and from 
Lang-Firsov strong coupling theory. Agreement of both per- 
turbative approaches with the exact numerical result is excel- 
lent. The MA result (not shown) is also in excellent agreement 
with the numerical results. (b)Effective mass m* /mso vs. A. 
Numerical results are compared with that from weak coupling 
perturbation theory, and agreement is excellent for low val- 
ues of A. Both exact and perturbative approaches show an 
enhanced effective mass with increasing spin orbit coupling. 



coupling A, with non-zero values of the spin orbit inter- 
action, Vs/t = 0.5 and Vg/t = 1.0; these are compared 
with the results from the Holstcin model with Vs/t = 0. 
Here the phonon frequency is set to be w^/t = 1.0, which 
is the typical value used in Ref.[ll|, and for each value 
of Vs, the ground state energy is compared to the cor- 
responding result for A = 0. The numerical results are 
compared with results from weak coupling perturbation 
theory and from Lang-Firsov strong coupling theory. 

In Figl^a), the ground state energy crosses over 
smoothly (at around A 0.8) from the dclocalizcd elec- 
tron regime to the small polaron regime. Note that there 
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FIG. 4. (color online) (a) Ground state energy Egs — Ed 
vs. A for Vs/t — 0,0.5, 1.0 and ujE/t = 0.1. Exact numerical 
results are compared with those from weak coupling perturba- 
tion theory (labeled Pert, in the fig) and Lang-Firsov strong 
coupling theory, (b) Ground state energy Egs — Eo vs. A 
in the weak and intermediate coupling regime, (c) Effective 
mass m* /mso vs. A. Numerical results are compared with 
those from weak coupling perturbation theory. 
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FIG. 5. (color online) Effective mass m* /mso vs. us/t for 
weak electron phonon coupling A = 0.064. In the inset the 
effective mass in the phonon frequency region near the adi- 
abatic limit is shown. It is clear that the effective mass is 
enhanced as spin orbit interaction decreases near the adia- 
batic limit. This is in agreement with the result inferred from 
the electron density of states shown in Fig. 2(a). 



is a slight dependence of the ground state energy on the 
spin orbit interaction. If we define AE — Egs — Eq, then 
AE(V s /t = 0.5) < AE{V s /t = 0) < AE(V s /t = 1.0) in 
the delocalized electron regime, which is in agreement 
with the weak coupling perturbation theory, though this 
is barely visible in the figure. In the small polaron regime, 
the ground state energy is shifted up by the spin or- 
bit interaction. This trend agrees with the results from 
Lang-Firsov strong coupling theory. For Vs/t = 0, the 
Lang-Firsov theory agrees very well with the numerical 
results, while as the spin orbit coupling Vs increases, the 
Lang-Firsov theory becomes less accurate for the same 
electron phonon coupling (e.g. if we look at A = 1.0, 
for Vs/t = 1.0, the difference between Lang-Firsov the- 
ory and exact numerical results is larger than that for 
Vs/t = 0). This is due to the fact that the bandwidth 
is increased by spin orbit interaction,— so the effective 
electron phonon coupling is decreased by spin orbit in- 
teraction. Better agreement with Lang-Firsov theory is 
achieved for larger values of A. In FigfJJb), the effective 
mass is enhanced by the spin orbit interaction in the de- 
localized electron regime, which is in agreement with the 
prediction from weak coupling perturbation theory. Here 
we have only shown results in the region Vs/t = ~ 1.0; 
for larger values of Vs/t the effective mass will be de- 
creased by the spin orbit interaction in the delocalized 
regimei^i In the small polaron regime, the effective mass 
will always be decreased by the spin orbit interaction. 
In Fig 21 we show the same results as Fig. [3] for a much 
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smaller phonon frequency u^/i = 0.1, which is closer to 
the adiabatic limit. In Fig|3ta), the ground state energy 
crosses over sharply (but still smoothly) from the delo- 
calized electron regime to the small polaron regime. If 
we use A c to describe the critical value for this sharp 
crossover, A c will be enhanced significantly by the spin 
orbit interaction. For Vs/t = 0.0, A c ~ 0.55, while for 
Vs/t = 5.0, A c ~ 1.55 from our numerical results. In 
Fig. QJb), in the delocalized electron regime the ground 
state energy is decreased by the spin orbit interaction 
AE{V s /t = 1.0) < AE(V s /t = 0.5) < AE(V s /t = 0.0), 
which is also in agreement with the weak coupling pertur- 
bation theory. For larger Vs/t the ground state energy 
will be increased^ in the delocalized electron regime. In 
the small polaron regime, the ground state energy will be 
increased by the spin orbit interaction, in agreement with 
the Lang-Firsov theory. In Figj4jc), the effective mass 
enhancement for different spin orbit interaction Vs/t is 
shown vs. electron phonon coupling strength, A. For 
Vs/t = there is a rather sharp crossover from the delo- 
calized electron regime to the small polaron regime^ 

Near the crossover point, the effective mass enhance- 
ment for the delocalized electron is around 1.4. For 
nonzero Vs/t < 1, near the crossover point, the effec- 
tive mass enhancement is higher, but still within the 
same order of magnitute as Vs/t = 0. This does not 
agree with the exact adiabatic limit w^/t = of the 
Rashba-Holstcin model which has been studied recently 
by Grimaldi^ based on a semiclassical method. He 
found that for nonzero spin orbit interaction Vs, the 
ground state will experience two phase transitions as the 
electron phonon coupling A is increased. The first tran- 
sition is from a delocalized electron to a large polaron, 
while the second one is from a large polaron to a small 
polaron. In FigfJJc), we observe only one sharp crossover 
from a delocalized electron to a small polaron. Our re- 
sults did not exclude the possibilities that a large polaron 
regime will be found for cjg/t < 0.1, although we find 
this possibility unlikely. A similar circumstance holds in 
the absence of a spin orbit coupling, where the adiabatic 
approximation gives rise to a single transition, while the 
quantum calculations results only in a crossover. Smaller 
values of uje can be explored, but quantum fluctuations 
become stronger for LUs/t < 0.1 and the problem is 
numerically expensive for intermediate electron phonon 
coupling. 

To obtain some insight for the polaron effective mass 
near the adiabatic limit, we resort to weak coupling per- 
turbation theory. In Fig|3] we observe an anomalous in- 
crease of the effective mass for small uje for nonzero Vs- 
However, the effective mass stops increasing as it reaches 
some finite number (around 1.2 and 1.1 for Vs/t = 0.5 
and 1.0, respectively), so this does not indicate a break- 
down of the perturbation theory. This result is confirmed 
by the MA results, as illustrated. This is also in agree- 
ment with results from the adiabatic limit. As shown in 
Fig. 2 of Ref.pl], for V s /t = 0.5 and 1.0 (his j/t = 1.0 
and 2.0), the electron is definitely in the delocalized elec- 
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FIG. 6. (color online) Effective mass m* /mso map as a func- 
tion of spin orbit interaction Vs/t and coupling constant A 
for LdE/t — 0.1 obtained with the momentum average approx- 
imation. 

tron regime for A = 0.064 (A = 0.1 in Rcf.pl]). Actually 
this anomalous increase of effective mass is caused by an 
increase in the value of the electron DOS at the bottom 
of the band, as shown in Fig. 2 and Eq. (|9|). Thus, 
for even smaller values of Vs/t, the anomalous mass en- 
hancement will increase further and perturbation theory 
will eventually break down. This is in agreement with 
the adiabatic limit results — as Fig. 2 of Ref.pl shows, 
for Vs/t ~ 0, the electron enters the large polaron regime 
for small A. As mentioned earlier, our results are consis- 
tent with crossovers rather than transitions. This can be 
also seen in Fig. [6] where we plot for completeness a map 
of the effective mass as a function of V s /t and A obtained 
by using the MA approximation for u^/t = 0.1. The ex- 
act results, while different in the details, show the same 
qualitative trends. 

In Fig. wc compare exact numerical results with 
both the momentum average method^ and with weak 
coupling perturbation theory, for different values of uie- 
In Fig.[7Ja), the ground state energy is shown as a func- 
tion of Vs/t, while in Fig. EJb), the effective mass is 
shown as a function of Vs/t. The MA method agrees 
well with the exact numerical results for w^/t = 1.0. 
For smaller values of oje (y>E/t = 0.1 and 0.2), the Mo- 
mentum Average approximation becomes less accurate 
and agrees more closely with weak coupling perturbation 
theory. This is similar to what happens for the Holstein 
model. Reasons for this quantitative failure of MA in the 
adiabatic limit are explained in Ref. Il2l . 



IV. SUMMARY 

In this paper we have studied the problem of a single 
electron coupled to oscillating ions, in the presence of a 
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approximate solutions must converge. Moreover, in many 
dilute semiconductor applications, the single electron re- 
sult may be the relevant regime required for understand- 
ing of the problem. 

The exact method of solution utilizes the Trugman 
method of solution^ through Lanczos diagonalization. 
The procedure for this is now well documented, and con- 
verges very quickly over a very wide parameter regime. 
The momentum average approximation^ also works very 
well over the entire parameter regime; there is a break- 
down for very low phonon frequencies. In this regime the 
adiabatic approximation^ provides a good qualitative 
picture. Weak coupling perturbation theory^ tends to 
be fairly accurate only for very small coupling strengths. 
Finally, strong coupling perturbation theory^ is very ac- 
curate in the small polaron regime. 

In weak coupling the presence of spin orbit coupling 
increases the effective mass of the electron coupled to 
Einstein phonons^ The effective mass is small to begin 
with, so in this regime the impact of spin orbit coupling is 
fairly minor. As the electron phonon coupling increases, 
and one enters the small polaron regime, the presence of 
spin orbit coupling has the opposite effect, as first noted 
with the momentum average approximation^ Since in 
this regime the effective masses can be quite large, spin 
orbit coupling can have a profound effect on the charac- 
teristics of the electron. 
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FIG. 7. (color online) (a) Ground state energy Egs — £is a 
function of spin orbit interaction Vs/t for uiE/t = 0.1, 0.2, 1.0 
at A = 0.32. (b) Effective mass m* /mso as a function of 
spin orbit interaction Vs/t for the same parameters. Exact 
numerical results are compared with those from momentum 
averaging methods and weak coupling perturbation theory. 
The Momentum Average approximation does not do as well 
for low phonon frequencies. 



spin-orbit interaction. This problem has become relevant 
for a variety of spintronics applications. Many previous 
treatments have addressed this problem with a finite den- 
sity of electrons, and have therefore necessarily required 
approximate theoretical methods for solution. The limit 
of only one electron, previously solved with weak cou- 
pling perturbation methods and with the momentum av- 
erage approximation, is amenable to exact solution as de- 
scribed here, and serves as a benchmark to which other, 



Appendix A: Density of states at the bottom of the 

band 



Expanding Eh,- around the minimum energy Eg, by 
defining k' x = k x ± arctan(-^), k' y = k y ± arctan(-^), 
we have 

£k,- — Eo 



0.5* 



[(l+V s 2 /t')(k'J+k' y ')±2k' x k' y ] 
(Al) 



To calculate the density of states at the bottom of the 
band, from the definition, we have 

D_(E Q + E\) = ——j [ dk x f dk y S(E + E 1 - e fe> _), 

J — 7T J ■ — TT 

(A2) 

where E\ is a small amount of energy above the bottom 
of the band, E . Around the four energy minimum points 
there are four small regions which will contribute to this 
integral. We choose one of them (and then times our 
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results by a factor of 4) and use the definitions of k' 
above instead of k, introduce a small cutoff k c which is 
the radius of a small circle around k m i n , thus the integral 
reads 



D_(^ + Bi)=4xl 



k'dk' 



(W 



xS{E 1 



0.5i 



[(l + V§/t 2 )+sm26]k' 2 ) 



f m 

„ [(l + V 2 /t 2 )+sin20} 



V2 1 

7T Vs' 



(A3) 



The derivation of the effective mass in the weak cou- 
pling approximation (Eq. ([20]) ) proceeds similarly. We 
begin with Eq. (fT8|) in the text for the self energy. For 
very small phonon frequency we need only focus on the 
lower Rashba band, s = — 1. Furthermore, the non- 
interacting electron energy can be expanded about a min- 
imum, as in Eq. (]Aip . Noting that there are four equal 
contributions coming from the four degenerate minima, 
we obtain 



£weak(w + id) 



■nXtujE 

i 



dk' / dkt 



i+- 



{i + {v s /m{k'i+k^)+2k' x k 



(A4) 



where a 2 = Eq+lue— u, and the integration is understood 
to be around a small disk located at one of the energy 
minima. Transforming to polar coordinates allows both 
the radial and angular integral to be done analytically; 
for the radial integral we keep only the dominant portion 
for small uje, and, after differentiation, we readily obtain 
the result quoted in the text (Eq. ([20]) ). 



Appendix B: Strong coupling limit 

To investigate the strong coupling limit using second 
order perturbation, we need to evaluate Eq. (|15[) . re- 
peated here for convenience: 



E 



(2 



N 

k- - E E 

ri TOT ^0.nl,n2,...=0,l,...oo <=i 

(rti,n 2 , —ntflph ® (c£ a \elT\^k,-)el ® |0) P J 



-TlTOT^E 



-t 2 e~ 2 9 

UJ E 



N 



E E 

nl,n2..=0 l =1 
™TOT^ n 



\A^\ 2 + \Atf 
nxoT 



(Bl) 



where A a is given a series of matrix elements (distinct for 
(7 =| and |). These turn out to give equal contributions, 



so we illustrate in some detail the result for A\ only. 
After some algebra, we obtain 



l^tl 2 = M-5)l 2 1 E w*G?)l s 

5— ±a:,±y 



(B2) 



where 



u e (±g) = (n.le^lO) 



I Jig}. 



(B3) 



and 



C+y 



D -j-ik x a 



(i 



V s 
— i 

"(1-1 



and 



sin (k y a) — i sin (k x a) 
J sin 2 (k x a) + sin 2 (fc y a) 



(B4) 



(B5) 



For each of the ug(±g) in Eq. (jB2|) it is to be understood 
that ng ^ 0, but all other ngi = for £' ^ £. Hence, in 
the 16 terms in Eq. (]B2[) . 12 will have all phonon num- 
bers equal to zero (other than ng); the other 4 will have 
both m and ng +x (or rig and ng_ x , etc.) not equal to zero 
in general. As already mentioned, the contribution from 
|A|J 2 is identical to that from |^| 2 , so this merely gives 
us a factor of 2 in Eq. (|Blj) . Moreover, translational 
invariance makes the contribution from each site iden- 
tical, so the sum over sites is trivially performed. This 
equation then becomes 



Efl = -^^{/(5 2 )(^) 2 +[/(2.9 2 )-/(5 2 )] [l+(^) 2 ]}, 



where 



(B6) 



= E ~—\ = Ei W - 7 - Inx 

n=l 

^e x /x[l + l/x + 2/x 2 + ...], (B7) 

and Ei{x) is the exponential integral and 7 « 0.5772 is 
Euler's constant. Eq. (]B6[) leads directly to Eq. ([TBI in 
the text. 
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